Extracting the qq condensate for light quarks 
beyond the chiral limit in models of QCD 
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It has recently been suggested [l[ that a reliable and unambiguous definition of the non- 
perturbative massive quark condensate could be provided by considering a non positive-definite 
class of solutions to the Schwinger Dyson Equation for the quark propagator. In this paper we show 
that this definition is incomplete without considering a third class of solutions. Indeed, studying 
these three classes reveals a degeneracy of possible condensate definitions leading to a whole range 
of values. However, we show that the physical condensate may in fact be extracted by simple fitting 
to the Operator Product Expansion, a procedure which is stabilised by considering the three classes 
of solution together. We find that for current quark masses in the range from zero to 25 MeV or so 
(defined at a scale of 2 GeV in the MS scheme), the dynamically generated condensate increases 
from the chiral limit in a wide range of phenomenologically successful models of the confining QCD 
interaction. Lastly, the role of a fourth class of noded solutions is briefly discussed. 

PACS numbers: 12.38.-t, 11.30.Rd, 12.38.Aw, 12.38.Lg 
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I. INTRODUCTION 



A remarkable aspect of strong coupling field theory is 
the possibility that masses can be largely, or even en- 
tirely, created by interactions. It is by this mechanism 
that a fermion mass gap is generated for the light flavours 
in QCD, leading to dynamical chiral symmetry break- 
ing (DCSB). Indeed, the mass gap generated far exceeds 
the scale of the current quark masses present in the La- 
grangian. The strong coupling induces long range qq cor- 
relations that polarise the vacuum. It is the scale of these 
condensates that determines the mass generated, and this 
gap persists even in the chiral limit. As long known this 
underpins much of QCD phenomenology. In this paper 
we investigate this mass generation using the Schwinger- 
Dyson equations. The aim is to extract the behaviour of 
the qq condensate beyond the limit of zero quark mass. 

The interest in the value of such a condensate arises 
in the context of QCD sum-rules. There the Operator 
Product Expansion (OPE) is used to approximate the 
short distance behaviour of QCD. In studying currents 
like that of ^7^(75)^, with qi = s and qj = u,d, the 
vacuum expectation values of uu, dd and ss operators 
naturally arise 0, H, H|- I n the chiral limit, the value of 
the qq condensate for the u and d quarks is well deter- 
mined to be —(235 ±15 MeV) 3 by experiment — in par- 
ticular from the low energy behaviour of 7T7t scattering || . 
However, in the OPE it is the value of the condensates 
away from the chiral limit that actually enters. Since 
the current masses of the u and d quarks are only a few 
MeV, the resulting condensate is expected to be close to 
its value in the chiral limit, but how close? For the first 
20 years of QCD sum-rules their accuracy was never suf- 
ficient for it to matter whether this difference was a few 
percent, 10% or even 20% effect. This equally applied to 
the estimate by Shifman, Vainshtein and Zakharov 



that the ss condensate was (0.8 ± 0.3) of the uu and dd 
values. It is the greater precision brought about by the 
studies of Refs. 0, H, [tj EH > for instance, that motivate 
the need to learn about how the qq condensate depends 
on the current quark mass. Indeed, Dominguez, Ram- 
lakan and Schilcher Q compute the ss condensate to be 
just (0.5 ± 0.1) times that for uu and dd. In the light 
of a better understanding [TJ E2, EH EH of strong cou- 
pling QCD how robust is this? First results were already 
presented in Ref. [lj|. 

For light quarks, u, d and s, studying the Schwinger- 
Dyson equation for the fermion propagator in the contin- 
uum is essential, until computation with large lattice vol- 
umes become feasible. Since the continuum Schwinger- 
Dyson equations can be solved for any value of the quark 
mass, they also provide a natural way to bridge the gap 
between lattice data at larger masses and the chiral limit 
of phenomcnological importance. Our primary focus is, 
of course, on QCD, but we shall draw on the NJL model 
where necessary. 



II. SCHWINGER- DYSON EQUATIONS 



Solving the Schwinger-Dyson Equations for QCD is an 
entirely non-trivial process, see Bili for re- 
views. They comprise an infinite tower of coupled inte- 
gral equations, that can be solved analytically only for 




FIG. 1: Schwinger-Dyson equation for the quark propagator 
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specific kinematical situations, for instance in the far in- 
frared and in the ultraviolet [l(| EE 0, EH • In general, 
the solutions have to be found numerically and then only 
after some form of truncation has been applied. To this 
end we require some suitable ansatz for the three-point 
functions in order to allow us to solve self-consistently 
for propagators. In the following sections we will use 
a number of truncation schemes that have been tested 
elsewhere [20l . |2~lI [22j], within which we may solve for 
the fermion propagator. Our chief aim is to calculate 
the mass function of the quark propagator for a range of 
current masses. Our starting point is the renormalizcd 
Schwinger-Dyson equation for the quark propagator as 
depicted in Fig. [TJ 



S^ip) = z 2 \s (0 \p) 
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Z 3 (27r)« 
x ltl S F (k)T v (k,p)D, lv {p-k) . (1) 

In the Landau gauge we can choose [23| Z\ = 1. The 
inverse propagator Sp l (p) is specified by two scalar func- 
tions A and M.: 



S~ F \v) = A(j> 2 ) (y + M{p 2 )) 



(2) 



While A is also a function of the renormalisation point fi 
and so strictly A(p 2 , /x 2 ), the quark mass function Ai(p 2 ) 
is renormalisation group invariant. Projecting out these 
two functions from Eq. (TTJ), we have two coupled equa- 
tions to solve. 

The operator product for the mass function has an 
expansion at large momenta given symbolically by: 



M{p 2 ) ~ m(p 2 ) + 



const 



m(p 2 )) 



(3) 



where the first term corresponds to the explicit mass in 
the Lagrangian, and the second to the lowest dimension 
vacuum condensate. For now we just show the momen- 
tum dependence given by the canonical dimensions, and 
leave for later the implications of the anomalous dimen- 
sions of QCD. Having computed the mass function using 
the Schwinger-Dyson equations, the essential problem is 
how to separate these two terms in Eq. with any ac- 
curacy if m is non-zero. We note that the mass function 
for a physically meaningful solution is expected to be 
positive definite. 

In the chiral limit, there exist three solutions for the 
mass function M(p 2 ). These correspond to the Wigner 
mode (the only solution accessible to perturbation the- 
ory), and two non-perturbative solutions of equal magni- 
tude generated by the dynamical breaking of chiral sym- 
metry. These we denote by: 



M( P 2 ) 



M w (p 2 ) = 
M±(p 2 ) = ±M°(p 2 ) 



(4) 



Such multiple solutions have been found in the context 
of QED4 by Hawes et al. 24]. One can ask whether 
analogous solutions exist in QCD as we move away from 
the chiral limit and what relevance they hold. Indeed, in 
a recent paper [l[ it was suggested that one could make an 
unambiguous definition of the massive quark condensate 
by taking a particular combination of these solutions. 
The existence of these is restricted to the domain: 



T> = {m : < m < m cr } , 



(5) 



where only the positive-definite solution exists beyond 
m cr . Chang et al. [l| found that inside the critical do- 
main, the solutions for both M + (p 2 ) and M~(p 2 ) exhibit 
the same running current-quark mass in the ultraviolet. 
In terms of Eq. (J3j) , this means both solutions have the 
same m term. Noting that the M~(p 2 ) solution had a 
condensate of opposite sign, they proposed a definition 
of the massive quark condensate given by: 



<7{m{n)) = lim ^(/i, A) N c tvo 



A— >oo 



rf 4 fc 

(2tt)4 



x - [S+ - S-] 



(6) 



where fj, denotes the renormalisation point. Since both 
propagators S± share the same asymptotic behaviour 
M.{p 2 ) — ► m(p 2 ) for large momenta, potential diver- 
gences in the integral cancel and the expression Eq. ^ 
is well defined. The resulting condensate cf(m(/x)) will be 
equal to the one for the physical M + solution provided 
S± have condensates of equal magnitude away from the 
chiral limit. However, we will show that this assumption 
is not correct. Morever, we find here that there is in fact 
an analogous solution Sw to the Wigner mode of the chi- 
ral limit. (This has also been noted in the revised version 
of This solution has the same ultraviolet behaviour 
of a running current quark mass, i.e. all three, M± and 
Mw, have a common m(p 2 ) in their OPE, see Eq. ([3]). 
Consequently, the combination 



S(0) = (2-P)S+- + 2(0 - 1) S 



(7) 



has its asymptotics controlled by the second term of the 
OPE, Eq. (2), for any f3. Thus, we can extend the def- 
inition of Chang et al. [1] to a family of condensates 
parametrised by [3: 

a(m( M ),/3) = lim Z 4 (^ A) N c tv D [ S(p) ,(8) 

where all SL, S+ and Sw are dependent upon the mo- 
mentum k 2 , the renormalisation point \i and the quark 
mass m(p). The choice (3=1 corresponds to the defini- 
tion of Eq. ©, with (3 = and (3 = 2 corresponding to 
two other natural choices. In fact (3 could take any value 
from —00 to +00 and so we can have a whole range of 
values for the massive quark condensate: all agreeing in 
the chiral limit. Consequently, the definition proposed 
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by Chang et at is far from unique, and does not provide 
a value for the condensate that corresponds to the phys- 
ical M + solution of interest. Eq.© merely defines the 
value for the difference of the condensates for M + and 
M~ . However, we will show that when combined with 
the OPE the 3 solutions will pick out a precise physical 
definition of this condensate. Before investigating this in 
the context of models of the QCD interaction, we shall 
draw analogy with the NJL model, within which a nat- 
ural definition of the massive quark condensate already 
exists. 



III. THE NAMBU-JONA-LASINIO MODEL 

Because of the complexity of QCD, it is often pru- 
dent to examine simpler systems exhibiting similar char- 
acteristics first. One such example is the Nambu-Jona- 
Lasinio (NJL) model. Though originally formulated to 
describe nucleon interactions in the pre-QCD era, the 
model can be reinterpreted by regarding the nucleons as 
quarks [25], HH- It can then be used to study bound 
states and so determine basic phenomenological quanti- 
ties of meson interactions. The NJL model shares the 
same symmetry structure as QCD, and is dominated by 
DCSB effects at low energies too. The Lagrangian for 
the NJL model with just two flavours of quarks with de- 
generate mass toq = m u = md is: 
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Cnjl = i)(x) (i Ip - TO ) ip(x) + C int 



(9) 



The interactions are given by a four-fermion contact 
term: 



OKxMx)) 2 + (*P(x)i l5 T a i;(x)) 2 } (10) 



with the two terms corresponding to the scalar and pseu- 
doscalar channels respectively. From this can be derived 
the so-called gap equation: 



to = too + i G^NcNf 



d 4 p 
(2-) 4 



tv D S(p) , (11) 



where the trace is over spinor indices in D-dimensions. 
At this point one has a choice of how to regulate the 
integrals. Since we are dealing with a non-renormalisable 
effective theory, our results will depend upon the cut-off 
used. One may introduce a non-covariant cut-off in the 
Euclidean 3-momentum, or employ a variety of covariant 
regularisation schemes such as a four-momentum cut-off, 
proper time or Pauli-Villars. The four-momentum cut- 
off is most closely related to the scheme we will employ in 
the later sections, so we choose this method. By inserting 
the form of the propagator into Eq. |[TT|) we arrive at 



too + raG^Nc 



d 4 k 



1 



(2tt) 4 k 2 - TO 2 



(12) 
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FIG. 2: The non-perturbative condensate as a function of the 
current quark mass for a selection of NJL couplings G^. 



where we have introduced a cut-off A on the 4- 
momentum. Rotating to Euclidean space we obtain 



to 



G-kNcITI 

2tt 2 



A 2 - to 2 log ( 1 



A2 

9 

771 z 



(13) 



The coupling G^ and cut-off, A, are fixed by fitting to 
experimental data. To do this we must calculate the 
order parameter associated with the breaking of chiral 
symmetry. Thus the chiral condensate is given by: 



/d A k 
-^-^ tx D S F {k 



(14) 



and for non-zero current masses we define: 
d 4 k 



(qq) 



-iN n 



(2tt) 



tr£> [S F {k;m) - S F (k;m )} 



(15) 

Working with explicitly massless quarks and employing 
a covariant cut-off in the Euclidean 4-momentum, the 
pion-decay constant is given by [27] : 



fl = N c m 2 







log 






\ TO / 



A 2 



A 2 



(16) 



By solving Eq. (JT3J) and Eq. (fTfj|) simultaneously and de- 
manding that f n = 93 MeV and - (qq) 1/3 = 235 MeV, 
we obtain A = 0.908 GeV and to = 265 MeV. These can 
be substituted into the mass gap equation of Eq. (fTB")) and 
thus we can solve for the coupling, finding G^ = 10.2. 

Our interest is to investigate the condensate's depen- 
dence on the current quark mass. In Fig. [2] we see that 
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FIG. 3: The quark condensate within the NJL model for the 
three solutions M ,w as functions of quark mass. 



the behaviour of the condensate depends upon the cho- 
sen coupling. Indeed, it either increases to a maximum 
then decreases, or is monotonically decreasing for larger 
couplings as the quark mass increases. Moreover, one 
finds that the condensate can become vanishingly small 
for sufficiently large masses (~ 600 MeV, not shown in 
Fig. HI). 

Now, within the NJL model we also have solutions cor- 
responding to M + , M~ and M w away from the chiral 
limit within some domain. The extent of the domain de- 
pends on the parameters of the model. With the favoured 
choice, m cr ~ 15 MeV in Eq. ([5]). One may then use 
Eq. (jTSJ) to calculate the massive quark condensate for 
these solutions individually, the result of which is shown 
in Fig. [3] Furthermore, using the definition, Eq. ©, gives 
a family of condensates, <f(m(/i) , [3) , for which the results 
are shown in Fig. |4] for (3 = {0,1,2}. What is clearly 
evident is that the condensates for the three mass func- 
tions M ± , M w are not equal. Drawing analogy with the 
Schwinger-Dyson solutions for QCD, one should therefore 
not expect the condensate of M ± to be equal in magni- 
tude as well as opposite in sign, as assumed by Chang et 
al. [l|. We show this to be the case in the next section. 



IV. PHENOMENOLOGICAL MODEL OF QCD 
INTERACTION 



We now turn to QCD. Rather than solving for the 
ghost and gluon system, one may employ some suitable 
ansatz for the coupling which has sufficient integrated 
strength in the infrared so as to achieve dynamical mass 
generation. There have been many suggestions in the 
literature [2(| [28| which have been extensively studied. 



Following the lead of Maris et al. [2Cj, |28| , we will employ 
an ansatz for g 2 D fl „(p — k) which has been shown to be 
consistent with studies of bound state mesons. We will 
consider other modellings in a later section. Since this 
simple model assumes a rainbow vertex truncation, the 
solutions are not multiplicatively renormalisable and so 
depend on the chosen renormalisation point. For com- 
parison with earlier studies we take this to be /i = 19 
GeV. We will scale results by one loop running to 2 GeV 
in the modified momentum subtraction scheme relevant 
to the Maris- Tandy model. Thus we use: 



l-^-D^^a^D^q) 
where the coupling is described by: 



a(q 2 ) = —Dq^M-q'/u; 2 ) 



(17) 



log fr+ (i + q yA 2 QCD 

x [l-exp(- g 2 /[4m?])] , 



(18) 



with 



m t = 0.5 GeV 

7ro = 12/(33- 2V ; ) 



A 



QCD 



= 0.234 GeV. 



Note that we work in the V/ = limit first since in 
Sect. [Vl]we will investigate the mass dependence of the 
condensate using a model derived from quenched lattice 
data [22I ]. The precise value of Aqcd is irrelevent for 
our current study, and we choose the parameter set to = 
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FIG. 4: The condensate defined by <r(/3) within the NJL 
model for values f3 = {0, 1,2} as a function of quark mass. 
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FIG. 5: Euclidean mass functions for different current masses, 
specified at /i = 1 9 GeV as labelled. The plot illustrates how 
on a log-log plot the behaviour dramatically changes between 
a current mass of and 3 MeV. These results are essentially 
the same as found by Maris and Roberts [2Sl |. 
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FIG. 6: Momentum dependence of the three solutions M (p 2 ) 
and M w (p 2 ) for a quark mass m(»=16 MeV, fi = 19 GeV. 



0.4 GeV, D = 0.933 GeV 2 in the range considered by 
Ref. 0. 

Solutions are obtained by solving the coupled system 
of fermion equations for A and M. of Eq. which we 
may write symbolically as: 



A{ P \ti) = Z 2 ((M,A)--£ D (p,A) , 
M(p 2 )A( P 2 ,fi) = Z 2 (M>A)Z m mn(M) + Ss(p,A) 



(19) 



The Eg and Srj correspond to the scalar and spinor pro- 
jections of the integral in Eq. (JTJ) . For massive quarks we 
obtain the solution M + by eliminating the renormalisa- 
tion factors Z 2 , Z m via: 



Z 3 {n,A) = l + H D (jt,&) , 
Z m (p,A) = 



1 



(/z,A) 



(20) 



Z 2 (n,A) Z 2 (p,A)mR(fi) 



The resulting momentum dependence for different values 
of rriR are shown in Fig. [5] Our purpose is to define the 
value of the qq condensate for each of these. 

At very large momenta the tail of the mass function 
is described by the operator product expansion Eq. ([3]). 
Crucially, for a given wlr we obtain the M~ and M w 
solutions by inserting the same Z 2 and Z m found for the 
M + solution. This ensures that differences in the dynam- 
ics of the three systems do not influence the ultraviolet 
running of the current-quark mass in the context of the 
subtractive renormalisation scheme used by Maris and 
Tandy. The iteration process is performed using New- 
ton's method. For the solution M w this is mandatory, 
since it corresponds to a local maximum of the effective 
action and is therefore not accessible using the conven- 
tional fixed point iteration scheme. A representative ex- 
ample of the solutions is shown in Fig. [5J 



The value of the critical mass is model-dependent, and 
is summarised in Table HI Chang et al. [1] imbue this 
critical mass with some significance for the dynamics of 
QCD. However, criticality does not feature in the physical 
solution M + , which exists for all values of m q . It only 
occurs in the M~ and M w solutions, which appear in 
a strongly model-dependent region. Consequently, we 
find little evidence of criticality being important to the 
mass generation in QCD. We will comment again on this 
when we consider more sophisticated vertex structures in 
a later section. 

The definition S((3) of Eq. ([7]) gives for each (3 a mass 
function, for which the first term in the OPE, Eq. Q 
vanishes, and so is controlled entirely by the condensate 
term. However, as with the analogous NJL model, we 
have an infinite set of ambiguous definitions of the quark 
condensate, one for each value of (3, each of which agrees 
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TABLE I: How the critical mass that defines the domain of so- 
lutions Eq. ((5)1 depends on the number of quark flavours, Nf , 
on the gluon range parameter ui, in the Maris- Tandy model. 
This critical mass is listed at two different renormalization 
scales, 19 GeV of Ref. [IJ and 2 GeV for ease of comparison 
with other works in a momentum subtraction scheme. 
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FIG. 7: Renormalisation point independent quark condensate 
as a function of m q as defined by Eq. ([8]) for three values of 
/3, showing how they are quite different despite the solutions 
having the same running current-mass and being equal in the 
chiral limit. 



with the chiral condensate in the limit m q — > 0. This 
ambiguity seen in Fig. [7] arises because, although each 
solution exhibits the same leading logarithmic behaviour 
in the ultraviolet limit, the condensates for each are not 
equal in magnitude, cf. Eq. |3]). Indeed, the solutions 
M~ . M w have negative condensates, but we cannot di- 
rectly use combinations of these mass functions to form 
a well-defined and unique condensate that coincides with 
the true condensate contained within M + (p 2 ). 



V. EXTRACTING THE CONDENSATE 

At very large momenta the tail of the mass function is 
described by the operator product expansion of Eq. ([3]). 
For QCD, let us introduce the appropriate anomalous 
dimension factors explicitly, so that 



M(p 2 ) asym = m [log (pVa?)] " 



2?r 2 lm C 
3 p 2 



log (p 2 /A 2 - 



(21) 



where m is related to the quantity mfl(/x) via some renor- 
malisation factors. This provides an excellent represen- 
tation of all our solutions. If we included the expression 
to all orders then the scales A 1; and A2 would both be 
equal to Aqcd- However, the leading order forms in 
Eq. (|21[) absorb different higher order contributions into 
the two terms and so Ai and A2 are in practice different, 
as we will discuss below. For large masses the conden- 
sate piece, C, is irrelevant and so it is the leading term 
that describes the mass function well. In contrast in the 



chiral limit, m = and so the second term of the OPE 
describes the behaviour of the mass function. This then 
accurately determines the scale A2. Indeed, its value is 
equal to Aqc d ■ We can then easily extract the renormal- 
isation point independent condensate, C = — (qq), from 
the asymptotics — see Fig. 5 for the chiral limit. 

In this latter case, strictly in the chiral limit, we may 
also extract the condensate via: 



- (qq) = Z 2 (/x, A) Z m (zx, A) N c tr D 



d 4 k 



S (k,n) 



(22) 

where (qq) ^ is the renormalisation dependent quark con- 
densate. At one-loop, this is related to the renormalisa- 
tion point independent quark condensate: 




(23) 



which we compare with the asymptotic extraction to 
good agreement. 

However, for small quark masses, where the conden- 
sate is believed to play a sizeable role, we cannot ap- 
ply Eq. (|2"2")) . since it acquires a quadratic divergence, cf. 
Eq. ([3]). Indeed, it is the elimination of this that inspired 
the original Eq. ([6]) and later Eq. ([8]) , which we have seen 
lead to a wholly ambiguous definition of the physical con- 
densate. Nevertheless, one can attempt to fit both terms 
of the OPE in Eq. (j2~Tj) to the tail of the mass function, 
M + for instance. While a value for the condensate can 
then be extracted, this procedure is not at all reliable 
because of the difficulty in resolving the two functions in 
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FIG. 8: Condensate extracted through simultaneous fitting of 
the three solutions to the fermion mass- function in the Maris- 
Tandy model with Nf = and ui — 0.4 GeV as functions 
of the current quark mass defined at 2 GeV in a modified 
momentum subtraction scheme. 
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FIG. 9: Condensate for Maris-Tandy Model with Nf = 4, 
Lu — 0.4 GeV as a function of current quark mass defined at 
2 GeV, to be compared with Fig. [8] 



the OPE from one another and in fixing the appropriate 
scales, Ai and A2. 

It is to this last point which we now turn. Instead of 
one single solution, we now have three solutions to the 
same model, each with identical running of the current- 
quark mass (the first term in Eq. I|21|) ) in the ultraviolet 
region and differing only by their values of the conden- 
sate. Thus it is possible to fit Eq. (f2"Tj) simultaneously to 
the three mass functions M ± , M w . The scales Ai and 
A 2 are determined separately for each value of the cur- 
rent mass. Remarkably, within the given model, they are 
exactly the same for the range of quark masses we con- 
sider with Ai ~ 2Aqcd and A2 ~ Aqcd respectively. 
The condensates and C w are then determined in an 
accurate and stable way. For this to work the solutions 
have to be found to an accuracy of 1 part in f 8 . This fit- 
ting is performed using a modified Levenberg-Marquardt 
algorithm with appropriate weights added to give better 
behaviour at large momentum in accord with perturba- 
tion theory. The results for the phenomenological model 
employed here are given in Fig. [8l The error bars re- 
flect the accuracy with which the mass functions, repre- 
sentable by two terms in the OPE expression, Eq. (f2"Tj) . 
are separable with the anomalous dimensions specified. 
In contrast to the condensate defined by Eq. © , we find 
that in the limited mass range investigated, the conden- 
sate increases as a function of m q . At the critical point 
m cr (ii = 2 GeV)=44 MeV, we find the ratio for the con- 
densate to the chiral limit with Nf — to be (Fig. [S]): 



m=50Mev/ m— 



= 1.24 



(24) 



To estimate the errors in this determination, we can form 
combinations of these condensates in the same way as de- 
fined in Eq. ffSJ) , favourably reproducing the same results 



of Fig. [71 In Fig. [9] is a similar plot with Nf = 4 and 
lu = 0.4 GeV, illustrating how m cr changes compared 
with Fig. M 

We see that within errors the condensate is found to 
increase with quark mass. This rise at small masses was 
anticipated by Novikov et al. [30| combining a pertur- 
bative chiral expansion with QCD sum-rule arguments. 
That the chiral logs relevant at very small m q are barely 
seen is due to the quenching of the gluon and the rain- 
bow approximation of Eq. ([T]). As we will show in the 
next section when we model more complex interactions, 
including matching with the lattice, this effect remains 
small. 



VI. MORE SOPHISTICATED MODELS OF QCD 
INTERACTION 



We now consider the consequences of using more so- 
phisticated vertex structure for the quark-gluon interac- 
tion in the quark Dyson-Schwinger equation, Fig.[TJ The 
first framework we study is a truncation scheme intro- 
duced in [U [3l| . It involves replacing the bare quark- 
gluon vertex of Sect. lIVI with the Curtis-Pennington (CP) 
vertex [32j, thus ensuring multiplicative renormalizabil- 
ity for the fermion propagator. In the Yang-Mills sector 
of QCD ansatze for ghost and gluon interactions have 
been introduced, which enable a self-consistent solution 
for the ghost and gluon propagators. The second scheme 
we shall investigate is an ansatz for the quark-gluon ver- 
tex, which has been fitted to lattice results, and was pre- 
viously employed in Ref [13, • 



Continuum studies: CP vertex 

In this truncation scheme we use explicit solutions for 
the Dyson-Schwinger equations for the ghost and gluon 
propagators, given diagrammatically in Fig. 1101 This sys- 
tem of equations has been solved numerically in [34| . The 
relevant ansatze for the ghost-gluon and trip le-g luon ver- 
tices have been discussed in the literature 0, [U [34| . 
The solutions for the ghost and gluon propagator 



D G (p) 



G(p 2 



p fl Pu\ Z(p 2 ) 



can be represented accurately by 



Z(P 2 ) 



G(P 2 ) 



, 2 \\ 1+25 

a(/i) 
a{p 2 ) 



R 2 (p 2 ), 
R-\p 2 ), 



(25) 
(26) 

(27) 
(28) 



with 



8 



R(p 2 ) = 



C ( P 2 /Af. M r+rf( P 2 /Af. M ) 2 - 



(29) 



with the scale Aym — 0.658 GeV, the coupling a(fi) = 
0.97 and the parameters c = 1.269 and d = 2.105 in 
the auxiliary function R(p 2 ). The quenched anomalous 
dimension 7 of the gluon is related to the anomalous di- 
mension 5 of the ghost by 7 = —1 — 25 and 5 — —9/44 for 
N f = 0. The infrared exponent k = (93 - \/1201)/98 « 
0.595 Il6j |. The running coupling a(p 2 ) is defined via the 
nonperturbative ghost-gluon vertex, 



a(p 2 ) = a(p)G 2 (p 2 )Z(p 2 ) 



(30) 



and can be represented by 
1 



1 - 
4tt 



p 2 /Km 
1 



a(0)+p 2 /A2, M x 
1 



/3o Vln(p 2 /A 2 



V 2 I^ YM 1 



(31) 



The value a(0) « 8.915/AT C is known from an analytical 
infrared analysis [l6|. 

In the quark-DSE Eq. fTJ) we use the solution of 
Eq. (|27p for the gluon propagator together with an ansatz 
of the form 

r v (q,k) = V„ Abel (p,q,k)W- Abel (p,q,k), (32) 

where p and q denote the quark momenta and k the gluon 
momentum. The ansatz factorises into an Abelian and 
an non-Abelian part which are specified and discussed 
in detail in [2l|. Here we only need to remark that the 
Abelian part V* bel is identical to the CP vertex [HJ. This 
construction carries further tensor structure in addition 



to the 7^-piece, which makes it an interesting ansatz in 
comparison with the simple model in Sect. ITVl The cor- 
responding numerical solutions for the quark propagator 



are discussed in detail in 21|. Here we are only inter- 



ested in the chiral condensate as a function of the current 
quark mass. The corresponding results can be found in 
Fig. QT] Despite the complicated tensor structure of the 
vertex Eq. Q32p we find similar results for the condensate 
as previously. We were able to extract the condensate 
from all three solutions and M w , again for a re- 
stricted region of m q < m cr . The physical condensate 
rises again slightly for small current quark masses and 
bends down for larger ones. The critical value of m cr is 
found to be 20 MeV at fj, = 19 GeV. This corresponds to 
30 MeV at /x = 2 GeV. With the parameters of the 
condensate in the chiral limit is 270 MeV, rather than 
the phenomenological 235 MeV we have used. However, 
its dependence with quark mass hardly depends on this 
exact value and so at the critical point we obtain the 
ratio 



(W>m= 



m=30 MeV 



/(m) m=0 = 1-175 



(33) 



As with the phenomenological model considered in 
Sect. II VI we find a considerable increase of the chiral con- 
densate with the current quark mass. 



Lattice Model 

The third model we investigate has been defined in 
[22| . The idea is to solve the coupled system of gluon, 
ghost and quark Dyson-Schwinger equations on a com- 
pact manifold with periodic boundary conditions, similar 
to lattice QCD. For the vertices in the Yang-Mills sector 
the same truncation scheme as in the last section is em- 
ployed. However, for the quark-gluon vertex an ansatz 






FIG. 10: Dyson-Schwinger equations for the gluon and ghost 
propagator. Filled circles denote dressed propagators and 
empty circles denote dressed vertex functions. 
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has been specified such that lattice results for the quark 
propagator have been reproduced on a similar manifold. 
Solving the system also in the infinite volume/continuum 
limit one can then study volume effects in the pattern of 
dynamical chiral symmetry breaking [22I [33| • 

In the infinite volume/continuum limit, i.e on R , the 
solutions for the ghost and gluon propagator are given by 
Eq. ||27]) and Eq. (J25J). The ansatz for the quark-gluon 
vertex is 

v v {k,n 2 ) = 7i ,r 1 (fc 2 )r 2 (fc 2 , M 2 )r 3 (A : 2 , / i 2 ) (34) 

with the components 



as a convenient infrared cutoff for the logarithms. The 
quark mass dependence of the vertex is parametrised by 

a{M) = 



1 + a 2 M(C 2 )/A Q 



CD 



a 3 M*(e)/Al CD 



(38) 

where M(£ 2 ) is determined during the iteration process 
at £ = 2.9 GeV. The parameters have been fitted to 
lattice results using a staggered [35| and an overlap action 
[36] ]. Here we use only the fit to the overlap quark; the 
corresponding parameters are given in table [Til Explicit 
solutions for the quark propagators are discussed in [13] . 

This modelling provides a link between lattice results 
and the continuum. The vertex function, as defined by 
Eqs. (35-37), which is currently fitted to the lattice data, 
involves several different scales: A 9 , A Y m and Aq C d- In 
the continuum this results in a failure to reproduce the 
perturbatively determined anomalous dimensions in the 
OPE, Eq. (f2"Tj) . The extraction of the tiny condensate 
term in the OPE is very sensitive to these anomalous 
dimensions. Combined with the effect the uncertainties 
in determining the multiple scales has on Ai and A 2 in 
Eq. (|2ip leads here to much larger errors than in our pre- 
viously modellings. Nevertheless, for the massive con- 
densate we again find solutions similar to the previous 
sections. The condensate corresponding to the solu- 
tions are given in Fig. [T2] and the critical value m cr is 22 
MeV, see Table IIII1 The large errors point to the need 
for further studies of matching lattice on a torus to the 
continuum, if we are to extract reliable infinite volume, 
continuum quantities like the quark condensate. That 
is for the future. The ratio of the condensates is here 

(3<?) ro=2 2Mev/ (??>m=0 = L17 ' This is a little lower 

than for the previous models of Eqs. (24,33). However, 
uncertainties in extraction are considerably larger. Nev- 
ertheless the ratio is still bigger than one. 



T 1 {k 2 ) 



ln(fc 2 /A 2 



QCD 



r 2 (fc',/i') = G(k 2 ,^) g(c 2 ,m 2 ) W) 

x h [ln(fc 2 /A 2 +r)] 1+5 



(35) 



(36) 



„ „ „ a(M) + k 2 /K 2 nrn 

r 3 (fc 2 , M 2 ) = z 2 {^ 2 ) \' 2/ ' QCD , (37) 

1 + k 2 /A^ QCD 

where <5 = —9/44 is the (quenched) one-loop anomalous 
dimension of the ghost, j m — 12/33 the corresponding 
anomalous dimension of the quark and r — e — 1 acts 





h 


(GeV) 


(GeV) 


0,1 


«2 


03 


overlap 


1.31 


1.50 


0.35 


25.58 


3.44 


2.23 



TABLE II: Parameters used in the vertex model, Eqs. (|34I37[1 . 



Model 


CP Lattice 


m cr (fi = 19 GeV) 
m cr (n = 2 GeV) 


20 16 [MeV] 
30 22 [MeV] 



TABLE III: The critical mass for our quenched CP and Lat- 
tice model. 



VII. CONDENSATE BEYOND THE CRITICAL 
MASS RANGE: NODED SOLUTIONS 

We see in Figs. [H U that the M~ and M w solutions 
bifurcate below m cr ~ 43.4(44.0) MeV with lo = 0.4 GeV 
for Nf = 0(4) respectively. But what about the value 
of the condensate for the physical solution M + beyond 
the region where M~ and M w exist, i.e. mji(ii) > m cr ? 
Having accurately determined the scales Ai and A 2 in the 
OPE of Eq. (p?Tj) in the region where all 3 solutions exist, 
we could just continue to use the same values in fitting 
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FIG. 13: Momentum dependence of the 4 solutions for the 
fermion mass- function in the Maris- Tandy model with m = 20 
MeV at n = 19 GeV, N f =4, lo = 0.4 GeV. 



the physical M + solution alone and find its condensate. 
Unfortunately, this would make it difficult to produce 
realistic errors as the quark mass increases. 

However, as soon as one allows for solutions for the 
fermion mass-function that are not positive definite, one 
exposes a whole series of variants on the solutions M~ , 
M w we have already considered. Thus there are noded 
solutions, which have also been discovered recently in the 
context of a simple Yukawa theory by Martin and Llancs- 
Estrada • These noded solutions are only accessible if 
sufficient numerical precision is used. For instance, repre- 
senting the dressing functions by a Chebyshev expansion 
does not provide the required accuracy with sufficient 
smoothness to reproduce the OPE form of Eq. (21). We 
illustrate this within the Maris- Tandy model, for instance 
with Nf = 4 and w = 0.4, in Fig. [T2] There the four so- 
lutions we have found are displayed. 

It is interesting to note that this noded solution is not 
limited to the same domain that restricts M~ and M w . 
These noded solutions do develop a singularity in M{p 2 ) 
beyond m = 51.4 MeV at /x = 2 GeV. However, this 
is compensated for by a zero in A(p 2 ), Eq. ©, until 
in = 66.3 MeV. Thus there exists a solution with a well- 
defined ultraviolet running of the quark mass exactly as 
the M + solution, as far as m = 66.3 MeV. While at small 
quark masses we have all four solutions, at larger masses 
there are still two. Consequently, we can confirm that 
the scales Ai and A2 of Eq. ([2~Tj) as fixed for m < m cr are 
still well-determined by our fit procedure for m > m cr 
and so deduce the condensates. Indeed, fitting the M + 
and M^ oded at each value of mfl(/z) with common scales 
in the OPE equation, Eq. (|21|1 allows the condensate for 
the physical solution to be found for much larger quark 
masses, as shown in Fig. 1141 Indeed, these fits confirm 



that Ai and A2 are independent of 771r(/x). We can then 
fit the remaining M + solutions shown in Fig. [5] to give 
the physical condensate shown in Fig. [14] for acceptable 
values of ui as determined by [29| . 

In Fig. [15] we scale the quark mass from fj, = 2 GeV in 
the (quark-gluon) MOM scheme by one loop running to 
the MS scheme at 2 GeV using the relationship between 
Amom and Ag? for 4 flavours deduced by Celmaster 
and Gonsalves [38[. In this latter scheme the strange 
quark mass is ~ 95 MeV as given in the PDG Tables [3!| • 
Within the range of the Maris- Tandy modelling of strong 
coupling QCD, we find the ratio of the condensates at the 
strange quark mass to the chiral limit is 



,1/3 

m(MS) =95 MeV 



/(qq) 



1/3 



1.1 ± 0.2) . (39) 



in a world with 4 independent flavours. Moreover, here 
all the quarks have the same mass and there is no mixing 
between different hidden flavour pairs. Elsewhere we will 
illustrate the change that occurs in solving the quark 
Schwinger-Dyson equations with 2 flavours of very small 
mass m u ^d and 1 flavour with variable mass. Of course, 
in the quenched case quark loops decouple and exactly 
replicate the results given here. 

What we have shown here is that there is a robust 
method of determining the value of the qq condensate 
beyond the chiral limit based on the Operator Product 
Expansion. Of course, as the quark mass increases the 
contribution of the condensate to the behaviour of the 
mass function, Fig. becomes relatively less important 
and so the errors on the extraction of the physical con- 
densate increases considerably. Nevertheless, the method 
is reliable up to and beyond the strange quark mass. Al- 
ternative definitions are not. 
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While the existence of multiple solutions to the fermion 

Schwinger-Dyson equation is essential for our method, 

only the M + solution has any physical significance and 

the others are mathematical curiosities. Those, like M~, 

M w and M^ ded , only exist in restricted domains. In 

contrast, the physical solution exists for all current quark 

masses, even if we cannot reliably extract the value of the 

corresponding condensate from the OPE. While it is clear 

that the radius of convergence of the chiral expansion of 

M + [40L l4lj in terms of the quark mass has a scale of 

1/3 

order of Aqcd or equally (—(qq)) , is this scale set by 
m cr of Eq. (5)? Chang et al. claim it is [1]. However, 
the bifurcation point for the unphysical solutions differs 
whether they are noded and not, cf. Figs. 9 and 14, and 
in turn each is highly model-dependent. This makes it 
difficult to claim that the value of m cr of Eq. (5) is the 
key parameter of the radius of the convergence of the 
chiral expansion for M + . 



VIII. SUMMARY 

Within the NJL model and the Schwinger-Dyson ap- 
proach to QCD, we have investigated the three incquiva- 



lent solutions, called M + , M~ and M , to the mass 
gap equation that exist within the interval T>(m) — 
{m : < to < m cr }. By ensuring each were solved using 
the same renormalisation conditions, we found that each 
solution exhibited the same running of the current-quark 
mass in the ultraviolet: differentiated solely by their in- 
frared behaviour and value of the quark condensate. 

Though it was not possible to define the condensate 
unambiguously by simply taking combinations of M + , 
M ~ and M w , the increased information available on the 
domain T>(ni) by having three solutions permits a reliable 
extraction of the condensate through simultaneous fitting 
of these to the OPE. In addition, we were able to obtain 
a fourth (noded) solution. This is only possible if the 
equations are solved to high numerical accuracy. Though 
this fourth solution violates the physical requirement of 
positivity, it has allowed us to extract the condensate 
beyond m cr and into the region of the physical strange 
quark mass. 

We have investigated a number of models for the strong 
coupling (infrared) behaviour of the quark-gluon inter- 
action of QCD and found in all cases that the con- 
densate, corresponding to the solution with a positive- 
definite mass function, increases moderately with current 
quark mass in the region under consideration. This is in 
contrast to the QCD sum-rule calculations of Refs. 6-8. 
Typically we find an increase of 30% from the chiral limit 
to a current MS mass of 100 MeV at a scale of 2 GeV. 
Only at still larger quark masses does the condensate 
significantly decrease. 
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